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ABSTRACT. It is shown that a constant magnetic field in 3+1 and 2+1 di- 
mensions is a strong catalyst of dynamical chiral symmetry breaking, leading to 
the generation of a fermion dynamical mass even at the weakest attractive inter- 
action between fermions. The essence of this effect is the dimensional reduction 
D ^ D — 2 in the dynamics of fermion pairing in a magnetic field. The effect 
is illustrated in the Nambu-Jona-Lasinio (NJL) model and QED. In the NJL 
model in a magnetic field, the low-energy effective action and the spectrum of 
long wavelength collective excitations are derived. In QED (in ladder and im- 
proved ladder approximations) the dynamical mass of fermions (energy gap in 
the fermion spectrum) is determined. Possible applications of this effect and its 
extension to inhomogeneous field configurations are discussed. 
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1. INTRODUCTION 



At present there are only a few firmly established non-perturbative phenomena 
in 2+1 and, especially, (3 + l)-dimensional field theories. In this paper, we will 
establish and describe one more such phenomenon: dynamical chiral symmetry 
breaking by a magnetic field. 

The problem of fermions in a constant magnetic field had been considered by 
Schwinger long ago In that classical work, while the interaction with the 
external magnetic field was considered in all orders in the coupling constant, 
quantum dynamics was treated perturbatively. There is no spontaneous chiral 
symmetry breaking in this approximation. In this paper we will reconsider this 
problem, treating quantum dynamics non-perturbatively. We will show that in 
3+1 and 2+1 dimensions, a constant magnetic field is a strong catalyst of dy- 
namical chiral symmetry breaking, leading to generating a fermion mass even at 
the weakest attractive interaction between fermions. We stress that this effect is 
universal, i.e. model independent. 

The essence of this effect is the dimensional reduction D ^ D — 2 in the dynamics 
of fermion pairing in a magnetic field: while at D = 2 + 1 the reduction is 
2 + 1— >-0 + l, atD = 3 + litis3 + l— ^1 + 1. The physical reason of this 
reduction is the fact that the motion of charged particles is restricted in directions 
perpendicular to the magnetic field. 

Since the case of 2+1 dimensions has been already considered in detail in Ref. 0, 
the emphasis in this paper will be on (3 + 1) -dimensional field theories. However, 
it will be instructive to compare the dynamics in 2+1 and 3+1 dimensions. 

As concrete models for the quantum dynamics, we consider the Nambu-Jona- 
Lasinio (NJL) model and QED. We will show that the dynamics of the lowest 
Landau level (LLL) plays the crucial role in catalyzing spontaneous chiral sym- 
metry breaking. Actually, we will see that the LLL plays here the role similar to 
that of the Fermi surface in the BCS theory of superconductivity [Q. 

As we shall show in this paper, the dimensional reduction D D — 2 is reflected 
in the structure of the equation describing the Nambu-Goldstone (NG) modes 
in a magnetic field. In Euclidean space, for weakly interacting fermions, it has 
the form of a two-dimensional (one-dimensional) Schrodinger equation at D = 



Here ^'(r) is expressed through the Bethe-Salpeter (BS) function of NG bosons. 



3 + 1 (D = 2 + 1): 




(1) 



A 



dxl dxl 
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(the magnetic field is in the +X3 direction, X4 = it) for D = 3 + 1, and A = 
X3 = it, for D = 2 + 1. The attractive potential V{r) is model dependent. In the 
NJL model (both at D = 2 + 1 and D = 3 + 1), V{r) is a (5-like potential. In 
(3 + 1) -dimensional ladder QED, the potential V{r) is 

where Ei{x) = — dtexp{—t)/t is the integral exponential function 0, « = |^ 
is the renormalized coupling constant and £ = \eB\~^^'^ is the magnetic length. 

Since — w^dyn plays the role of energy E in this equation and V{r) is an attractive 
potential, the problem is reduced to finding the spectrum of bound states (with 
E = — 'm-dyn < 0) of the Schrodinger equation with such a potential. More 
precisely, since only the largest possible value of m'^y^ defines the stable vacuum 
0, we need to find the lowest eigenvalue of E. For this purpose, we can use results 
proved in the literature for the one-dimensional {d = 1) and two-dimensional 
{d = 2) Schrodinger equation |^. These results ensure that there is at least one 
bound state for an attractive potential for d = 1 and d = 2. The energy of the 
lowest level E has the form: 

E{X) = -ml^{X) = -\eB\f{X) , (3) 

where A is a coupling constant (A is A = G in the NJL model and A = a in 
QED). While for d = 1, /(A) is an analytic function of A at A = 0, for = 2, it 
is non-analytic at A = 0. Actually we find that, as G —* 0, 

Kn=\^B\ ^^J^ , (4) 

where Nc is the number of fermion colors, in (2 + l)-dimensional NJL model [^, 
and 

'^''y^- n """Pi \eB\N,G ) ' 

where g = N^GA"^ /{Att"^), in (3+l)-dimensional NJL model. In (3+l)-dimensional 
ladder QED, mdyn is 



"^dyn = C\/\eB\ exp 



2 \2a 



where the constant G is of order one and a is the renormalized coupling constant. 
As we will show below, this expression for mdyn in QED is gauge invariant. 

It is important that, as we shall show below, an infrared dynamics, with a weak 
QED coupling, is responsible for spontaneous symmetry breaking in QED in a 
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magnetic field. This suggests that the ladder approximation can be reliable in 
this problem. However, because of the (1 + l)-dimensional form of the fermion 
propagator in the infrared region, there may be also relevant higher order con- 
tributions. As we shall show in this paper, there are indeed relevant one-loop 
contributions in the photon propagator. Taking into account these contributions 
(that corresponds to the so called improved ladder approximation), we get the 
expression for rridyn of the form with a a/2. We shall discuss the physics 
underlying this effect and the general status of the validity of the expansion in a 
in the infrared region in QED in a magnetic field in Sec. 9. 

As we will discuss in this paper, there may exist interesting applications of the 
phenomenon of chiral symmetry breaking by a magnetic field: in planar condensed 
matter systems, in cosmology, in the interpretation of the heavy-ion scattering 
experiments, and for understanding of the structure of the QCD vacuum. We will 
also discuss an extension of these results to inhomogeneous field configurations. 

The paper is organized as follows. In Section 2 we consider the problem of a 
free relativistic fermion in a magnetic field in 3+1 and 2+1 dimensions. We 
show that the roots of the fact that a magnetic field is a strong catalyst of 
dynamical chiral symmetry breaking are actually in this problem. In Section 3 
we show that the dimensional reduction 3 + 1 — 1 + 1 (2 + 1— s>0 + l) in the 
dynamics of fermion pairing in a magnetic field is consistent with spontaneous 
chiral symmetry breaking and does not contradict the Mermin-Wagner-Coleman 
theorem forbidding the spontaneous breakdown of continuous symmetries in less 
than 2+1 dimensions. In Sections 4-8 we study the NJL model in a magnetic field 
in 3+1 dimensions. We derive the low-energy effective action and determine the 
spectrum of long wavelength collective excitations in this model. We also compare 
these results with those in the (2 + l)-dimensional NJL model In Section 
9 we study dynamical chiral symmetry breaking in QED in a magnetic field. 
In Section 10 we summarize the main results of the paper and discuss possible 
applications of this effect. In Appendix A some useful formulas and relations are 
derived. In Appendix B the reliability of the 1/A^c expansion in the NJL model 
in a magnetic field is discussed. In Appendix C we analyze the Bethe-Salpeter 
equation for Nambu-Goldstone bosons in QED in a magnetic field. 
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2. FERMIONS IN A CONSTANT MAGNETIC FIELD 



In this section we will discuss the problem of relativistic fermions in a magnetic 
field in 3+1 dimensions and compare it with the same problem in 2+1 dimensions. 
We will show that the roots of the fact that a magnetic field is a strong catalyst 
of chiral symmetry breaking are actually in this dynamics. 

The Lagrangian density in the problem of a relativistic fermion in a constant 
magnetic field B takes the form 

£ = {zYD^ - m)^] , /X = 0, 1, 2, 3, (7) 

where the covariant derivative is 

D, = d,- zeA-* . (8) 
We will use two gauges in this paper: 

A^^* = -5,^Bx2 (9) 

(the Landau gauge) and 

= -l^i^iBx2 + ^S^2Bx, (10) 

(the symmetric gauge). The magnetic field is in the +0:3 direction. 
The energy spectrum of fermions is p|: 



En{k3) = ±^Jm^ + 2\eB\n + kl , n = 0,l,2, ... (11) 

(the Landau levels). Each Landau level is degenerate: at each value of the mo- 
mentum ks, the number of states is 



dNn = S,oL 



\eB\ dk'i 



1 9 

^ 27r 27r 

at n=0, and 

\eB\dk:, 
dN„. = S12L3- 



71 2% 



at n > 1 (where L3 is the size in the xs-direction and 5*12 is the square in the 
xia;2-plane). In the Landau gauge (||), the degeneracy is connected with the 
momentum ki which at the same time coincides with the X2 coordinate of the 
center of a fermion orbit in the a:iX2-plane P). In the symmetric gauge ([T0|), the 



5 



degeneracy is connected with the angular momentum J12 in the xiX2-plane (for 
a recent review see Ref. 0). 

As the fermion mass m goes to zero, there is no energy gap between the vacuum 
and the lowest Landau level (LLL) with n = 0. The density of the number of 
states of fermions on the energy surface with Eq = is 



yo = V 



£^0=0 



\eB\ 
4^ ' 



(12) 



£0=0 dEo 

where Eq = \ks\ and dNo = V^-^^ (here V = S12L3 is the volume of the 
system). We will see that the dynamics of the LLL plays the crucial role in 
catalyzing spontaneous chiral symmetry breaking. In particular, the density vq 
plays the same role here as the density of states on the Fermi surface in the 
theory of superconductivity 0]. 

The important point is that the dynamics of the LLL is essentially (l+l)-dimen- 
sional. In order to see this, let us consider the fermion propagator in a magnetic 
field. It was calculated by Schwinger |jl| and has the following form both in the 
gauge (^ and the gauge (pUj): 



S{x,y) = exp 



te 



(x-yrA^^\x + y) 



S{x - y) 



(13) 



where the Fourier transform of S is 



~S{k) 



(is exp 



,.2 ,2 ,2 tan(e5s) 



eBs 



m 



■[(A;V - 7 +m){l + Yl tan(e5s)) - kx7^(l + tan^(e5s))] (14) 



(here k_|_ = (ki,k2), 7_|_ = (71,72))- Then by using the identity itan(x) 
1 — 2 exp(— 2za;)/[l + exp(— 2za;)] and the relation |^ 



2)-(°+^)exp 



xz 



n=0 



(15) 



where L^{x) are the generalized Laguerre polynomials, the propagator S{k) can 
be decomposed over the Landau poles [llO| : 

k2 \ Dn{eB,k) 



S{k) = i exp 



\eB\ 



E(-i) 



n=0 



fin rv^ 



2\eB\n 



(16) 



with 



Dn{eB,k) 



{k^'Y - 7 + m) 



1 - z7Vsign(eS))LJ 2 



\eB\ 



[1 + z7S'sign(ei?))L„_i (27^)] + ^k'V + A;^7')^n-i (2- 



eB\ 



(17) 
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where L„ = and Lti = by definition. The LLL pole is 

~S^^\k) = ^ -p(-^)^^^f|^(l - ^lW^gn{eB)) . (18) 

This equation clearly demonstrates the (1 + 1) -dimensional character of the LLL 
dynamics in the infrared region, with -C \eB\. Since at m^, /cq, ^ 
\eB\ the LLL pole dominates in the fermion propagator, one concludes that the 
dimensional reduction 3 + 1 — > 1 + 1 takes place for the infrared dynamics in 
a strong (with \eB\ ^ m^) magnetic field. It is clear that such a dimensional 
reduction reflects the fact that the motion of charged particles is restricted in 
directions perpendicular to the magnetic field. 

The LLL dominance can, in particular, be seen in the calculation of the chiral 
condensate: 

(0|^^/^|0) = - limtr Six,y) = -j^ti [ d^k Ssik) 



4m 



m^ + kl + kl + kl'''''^^'^'^^ 



poo 

(27r)4 J A/A2 ^ . 
-\eB\^ r ^e-'-\otH\eBs\)^-\eB\^J\n^ + 0{m')), (19) 



where Ssik) is the image of S{k) in Euclidean space and A is an ultraviolet 
cutoff. As it is clear from Eqs. ([T6|) and (p!8|), the logarithmic singularity in the 
condensate appears due to the LLL dynamics. 

The above consideration suggests that there is a universal mechanism for enhanc- 
ing the generation of fermion masses by a strong magnetic field in 3+1 dimen- 
sions: the fermion pairing takes place essentially for fermions at the LLL and this 
pairing dynamics is (1 + 1) -dimensional (and therefore strong) in the infrared re- 
gion. This in turn suggests that in a magnetic field, spontaneous chiral symmetry 
breaking takes place even at the weakest attractive interaction between fermions 
in 3+1 dimensions. In this paper we will show the existence of this effect in the 
NJL model and QED. 

In conclusion, let us compare the dynamics in a magnetic field in 3+1 dimen- 
sions with that in 2+1 dimensions B. In 2+1 dimensions, we consider the four- 



component fermions |jTT|, connected with a four-dimensional (reducible) repre- 



sentation of Dirac's matrices 

-."-(? -''-(7 -D' -'^-(T -D- p«) 

The Lagrangian density is: 

C = \[^, (tYD^-m)^] , (21) 
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where = d^, - ieA"^^ with A"^^ = (0,-5x2,0) or A;f* = (0, -f xa, f xi). 
At m = 0, this Lagrangian density is invariant under the flavor (chiral) U{2) 
transformations with the generators 

To = /, Ti = 75, T2 = -7^ T3 = 7^7' (22) 

t 

where 7^ = i7''7-'^7^7^ — ( ^ t n V "^^^ mass term breaks this symmetry down 



-/ 0. 

to the U{1) X U{1) with the generators Tq and T3. The fermion propagator in 
2+1-dimensions is 0: 

5(x,i/)=exp(|(x-y)M-*(x + i/))^(x-?/) , (23) 

where the Fourier transform S{k) of S{x) is 



roo 

S(k) = / dsexp 

J1/A2 



2 7 2 2tan(ei?s) 

—zsm + iskr. — zsK 

ens 



■[(^7^ + m + (A;%^ - /c^^) tan(e5s))(l + 7^^ tan(e5s))] . (24) 
The decomposition S{k) over the Landau level poles now takes the following form: 

^(^.^^e.p(-Jj)|(-:r,^J:^ 

with 

k2 



Dn{eB,k) = (^7" + ^ 



l-27Vsign(e5))L„(2-^^ 



-(l + ^7Vsign(eB))L„_l 2 



k" 



2 



Ie5| 



4(fcV + A:V)TLi(2^) . (26) 



k' 



2 



Then Eq. (|25| ) implies that as m — 0, the condensate (0|'?/''?/'|0) remains non-zero 
due to the LLL: 

'm r ,n, cxp ( — k^ / 1 6,8 1 ) IcBl , 

(for concreteness, we consider m > 0). The appearance of the condensate in 
the flavor (chiral) limit, m — > 0, signals the spontaneous breakdown of f/(2) 
to t/(l) X f/(l) 0. As we will discuss in Section 5, this in turn provides the 
analyticity of the dynamical mass mdyn as a function of the coupling constant G 
at G = in the (2 + 1) -dimensional NJL model. 
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3. IS THE DIMENSIONAL REDUCTION 3+1 ^ 1+1 (2+1 ^ 0+1) 
CONSISTENT WITH SPONTANEOUS CHIRAL 
SYMMETRY BREAKING? 

In this section we consider the question whether the dimensional reduction 3+1 
1+1 (2+1 0+1) in the dynamics of the fermion pairing in a magnetic field 
is consistent with spontaneous chiral symmetry breaking. This question occurs 



naturaly since, due to the Mermin- Wagner- Coleman (MWC) theorem |jT2[, there 
cannot be spontaneous breakdown of continuous symmetries at D = 1 + 1 and 
D = + 1. The MWC theorem is based on the fact that gap less Nambu-Goldstone 
(NG) bosons cannot exist in dimensions less than 2+1. This is in particular 
reflected in that the (1 + 1) -dimensional propagator of would be NG bosons 
would lead to infrared divergences in perturbation theory (as indeed happens 
in the 1/Nc expansion in the (1 + 1) -dimensional Gross-Neveu model with a 
continuous symmetry [0). 

However, the MWC theorem is not applicable to the present problem. The cen- 
tral point is that the condensate (0| |0) and the NG modes are neutral in this 
problem and the dimensional reduction in a magnetic field does not affect the 
dynamics of the center of mass of neutral excitations. Indeed, the dimensional 
reduction D ^ D — 2 in the fermion propagator, in the infrared region, reflects the 
fact that the motion of charged particles is restricted in the directions perpen- 
dicular to the magnetic fleld. Since there is no such restriction for the motion of 
the center of mass of neutral excitations, their propagators have D-dimensional 
form in the infrared region (since the structure of excitations is irrelevant at 
long distances, this is correct both for elementary and composite neutral exci- 
tations). Q This fact will be shown for neutral bound states in the NJL model 
in a magnetic fleld, in the 1/Nc expansion, in Section 6 and Appendix B. Since, 
besides that, the propagator of massive fermions is, though (Z) — 2)-dimensional, 
nonsingular at small momenta, the infrared dynamics is soft in a magnetic fleld, 
and spontaneous chiral symmetry breaking is not washed out by the interactions, 
as happens, for example, in the (1 + 1) -dimensional Gross-Neveu model []13 . 



This point is intimately connected with the status of the space-translation sym- 
metry in a constant magnetic fleld. In the gauge (^, the translation symmetry 
along the X2 direction is broken; in the gauge (pUj), the translation symmetry 
along both the Xi and X2 directions is broken. However, for neutral states, all 
the components of the momentum of their center of mass are conserved quantum 



numbers (this property is gauge invariant) |[T^. In order to show this, let us in- 



^The Lorentz invariance is broken by a magnetic field in this problem. By the Z3-dimcnsional 
form, we understand that the denominator of the propagator depends on energy and all the 
components of momentum. That is, for D = 3 + 1, V{P) - {P^ - C^P^ - C'sP^y^ with 
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troduce the following operators (generators of the group of magnetic translations) 
describing the translations in first quantized theory: 



P.i = -7^ , P., = -7^ + QBx, , ^3 = (28) 
z oxi 1 0x2 ^ 0x3 

in gauge (||), and 

4^ = 1-^ - ^Bx2 , P^, = + -Bxi , P^, = (29) 

in gauge ( p!OD (Q is the charge operator). One can easily check that these operators 
commute with the Hamiltonian of the Dirac equation in a constant magnetic field. 
Also, we get: 

[P,„ P,,] = -QB , [P.,, P.3] = [P.„ P.3] = . (30) 

Therefore all the commutators equal zero for neutral states, and the momentum 
P = (Pi,P2,P3) can be used to describe the dynamics of the center of mass of 
neutral states. 
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4. THE NAMBU-JONA-LASINIO MODEL IN A MAGNETIC FIELD: 

GENERAL CONSIDERATION 



In this and the next four sections we shall consider the NJL model in a magnetic 
field. This model gives a clear illustration of the general fact that a constant 
magnetic field is a strong catalyst for the generation of a fermion dynamical 
mass. 

Let us consider the (3 + l)-dimensional NJL model with the x Uji{l) chiral 

symmetry: 

C = ^[ip, {t^DM + + i^tl'm , (31) 

where -D^ is the covariant derivative (§) and fermion fields carry an additional 
"color" index a = 1,2, . . . , N^. The theory is equivalent to the theory with the 
Lagrangian density 

C = l[iP, {tYDM - ^(a + Z7'vr)z^ - ^{^' + ■ (32) 

The Euler-Lagrange equations for the auxiliary fields a and tt take the form of 
constraints: 

a = -Giijip) , 7T = -Giipi-f^ifj) , (33) 

and the Lagrangian density (p2D reproduces Eq. (|3T|) upon application of the 
constraints (|33D . 

The effective action for the composite fields is expressed through the path integral 
over fermions: 

r(a, n) = f (a, ^) - ^ / dH^' + vr^) , (34) 



exp(ir) 



[diP][d,p] exp\ - I d'^xl^, [i^D^ -{a + z7^7r)}V^] 



exp < Tr Ln 



(cr + 27 tt) 



(35) 



I.e. 



f(cT, vr) = -zTr lj\i[iYD^ - {a + 27^71)] 



(36) 



As Nc ^ 00, the path integral over the composite fields a and vr is dominated 
by the stationary points of the action: 5T/5a = 5T/5tx = 0. We will analyze 
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the dynamics in this hmit by using the expansion of the action T in powers of 
derivatives of the composite fields. 

Is the 1/Nc expansion rehab le in this problem? The answer to this question is 
"yes". It is connected with the fact, already discussed in the previous section, 
that the dimensional reduction 3+1 1+1 by a magnetic field does not affect the 
dynamics of the center of mass of the NG bosons. If the reduction affected it, the 
1/Nc perturbative expansion would be unreliable. In particular, the contribution 
of the NG modes in the gap equation, in next-to-leading order in 1/Nc, would lead 
to infrared divergences (as happens in the (l + l)-dimensional Gross-Neveu model 
with a continuous chiral symmetry [0). This is not the case here. Actually, 
as we will show in Appendix B, the next-to-leading order in l/N^ yields small 
corrections to the whole dynamics at sufficiently large values of N^. 
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5. THE NJL MODEL IN A MAGNETIC FIELD: 
THE EFFECTIVE POTENTIAL 



We begin the calculation of the effective action T by calculating the effective 
potential V. Since V depends only on the Ui^{l) x t/R(l)-invariant = + tt^, 
it is sufficient to consider a configuration with tt = and a independent of x. So 
now r((T) in Eq. (^) is: 

T = -i TrLn{iD - a) = -i In Det(iD - a) , (37) 

where D = t'^-D^. Since 

Det(iZ) -o-) = Det(7^(iZ) - a)-f^) = Det{-iD - a) , (38) 

we find that 

f(a) = -^Tr[Ln(zi) - a) + Ln(-iZ) - a)] = -^Tr Ln(L)2 + a^) . (39) 

Therefore r(cr) can be expressed through the following integral over the proper 
time s: 

°° ds 



Y{a) = --Tt Ln{b^ + a^) = f d^x ( —ii{x\e-''^^'+^'\x) (40) 
2 2 J Jo s 



where 



= D,D>^ - -YYF;:' = D.D'^ + le-fWB . (41) 



The matrix element {x\e^^^^'^^"'^'^\y) was calculated by Schwinger |]l|: 

^^^^-is(b^+a^) ^ ^-isa^ ^^|g-isD,.DM [cos(e5s) + . sin(eSs)] 



-i,,, [g^^ cot(e5s) + 7^7'e5s] , (42) 



where 



{A-KS 



Sci = ef ATdz>^-Ux-yUg^^ + 

Jy As 



2 \u^i 
cxt/ 



■[l-eBscot{eBs)]){x-y)^ . (43) 

Here the integral Af^^dz'^ is taken along a straight line. Substituting expression 
(0) into Eq. (|g), we find 



Nr r r°° ds 



f((T) = ±4- d^x 4e"""'e5cot(eSs) . (44) 
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Therefore the effective potential is 



~, Nr ds 



V{p) = ^ + V{p) = ^ + 



— e"'^ e5coth(e5s) 



2G ' ' '''' 2G 87r2 7i/a2 
where now the uhraviolet cutoff has been exphcitly introduced. 



(45) 



By using the integral representation for the generalized Riemann zeta function ( 

i, 



(46) 



which is valid at /i > 1, and analytically continuing this representation to negative 
values of p, we can rewrite Eq. ( ^5|) as 



V{p) 



p" 



3^4 



2G 87r2 

(pA)2 + ^ \n{Mf + ^(1 - 7 - In 2) + ^ In 



Ac' 



+ 1 



0(1/A) 



(47) 



where the magnetic length £ = \eB\ ^1"^, C^'[—\,x) = -^Ci^^^)\u=~v 7 
0.577 is the Euler constant. [| 



The gap equation dV/ dp = is 



pA^--l 



-p^ In ^ — ^ + + In 

2 47r 



+ 2rVinr 



pH^ 



0(1/A) 



(4J 



where the dimensionless coupling constant g = N^GA"^ /{Att"^). In the derivation 
of this equation, we used the relations P] : 



d_ 

dx 



((z/,x) = -uC{u + l,x) 



(49) 



■^ai',x)l=o= Inr(x) - -ln27r , ({0,x) = --x. (50) 

^In this paper, for simplicity, we consider the case of a large ultraviolet cutoff: fP, 
where p is a minimum of the potential V. Notice, however, that in some cases (as in the 
application of the NJL model to hadron dynamics |^), and may be comparable in 
magnitude. 
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As i? ^ ^ oo), we recover the known gap equation in the NJL model |P, 



pA^^--lj= Vln^. (51) 

This equation admits a nontrivial solution only if g is supercritical, g > Qc = 1 
(as Eq. (|32|) implies, a solution to Eq. (^), p = a, coincides with the fermion 
dynamical mass, a = m^y^, and the dispersion relation for fermions is Eq. (pi]) 
with m replaced by a) . We will show that the magntic field changes the situation 
dramatically: at 7^ 0, a nontrivial solution exists for all g > 0. 

We shall first consider the case of subcritical g, g < gc = 1, which in turn can be 
divided into two subcases: a) g <^ gc and h) g ^ gc — Q (nearcritical g). Since 
aX g < gc = 1, the left-hand side in Eq. (|^) is positive and the first term on the 
right-hand side in this equation is negative, we conclude that a nontrivial solution 
to this equation may exist only at p^ln(A£)^ ^ £~^ln(p£)~^. Then we find the 
solution at (7 -C 1: 

One can check that the dynamical mass is mostly generated in the infrared re- 
gion, with k < = \eB\^^'^, where the contribution of the LLL dominates. In 
order to show this, one should check that essentially the same result for mdyn 
is obtained if the full fermion propagator in the gap equation is replaced by the 



LLL contribution and the cutoff A is replaced by J\eB\ (see Section 8 below). 



It is instructive to compare the relation (|5^) with the relations for the dynamical 
mass in the (1 + 1) -dimensional Gross-Neveu model [0, in the BCS theory of 



superconductivity [^], and in the (2 + 1) -dimensional NJL model in a magnetic 
field §. 

The relation for m^y^^ in the Gross-Neveu model is 



2 



2n \ 



^5yn = A exp(^-^^^j (53) 

where C-"^ is the bare coupling, which is dimensionless at D = 1 + 1. The 
similarity between relations (^) and (^) is evident: \eB\ and Ic-BIG in Eq. (^) 
play the role of an ultraviolet cutoff and the dimensionless coupling constant in 
Eq. (p3D, respectively. This of course reflects the point that the dynamics of 
the fermion pairing in the (3 + l)-dimensional NJL model in a magnetic field is 
essentially (l + l)-dimensional. We shall return to the discussion of the connection 
between this dynamics and that in the Gross-Neveu model in Sec. 7. 

We recall that, because of the Fermi surface, the dynamics of the electron in 
superconductivity is also (1 + l)-dimensional. This analogy is rather deep. In 
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particular, the expression ( p2D for mdyn can be rewritten in a form similar to 
that for the energy gap A in the BCS theory: while A ~ uoo exp(— const. /z/^^G^), 
where ujd'is the Debye frequency, Gs is a coupling constant and up is the density 



of states on the Fermi surface, the mass mdyn is mdyn ~ y exp(— l/2G'z/o), 
where the density of states z/q on the energy surface ii^ = of the LLL is now 
given by expression (0) multiplied by the factor Nc. Thus the energy surface 
E = Q plays here the role of the Fermi surface. 



Let us now compare the relation ( |52D with that in the (2 + 1) -dimensional NJL 
model in a magnetic field, in a weak coupling regime. It is 0: 

<~ = I'BI^^ (54) 

While the expression (^) for m^^^ has an essential singularity at G = 0, m-dyn in 
the (2 + l)-dimensional NJL model is analytic at G = 0. The latter is connected 
with the fact that in 2+1 dimensions the condensate (OlV'V'IO) is non-zero even 
for free fermions in a magnetic field (see Eq. (0)). Indeed, Eq. ( |33D implies 
that mdyn = (0|o-|0) = -G(0|V^^|0). From this fact, and Eq. (g^), we get 
the relation (^4|), to leading order in G. Therefore the dynamical mass mdyn is 
essentially perturbative in G in this case. As we will see in Section 8, this is 
in turn intimately connected with the fact that, for D=2+l, the dynamics of 
fermion pairing in a magntic field is (0 + 1) -dimensional. 

Let us now consider near-critical values of with A^(l — g)/gp^ ^ In(Ap)^. 
Looking for a solution to the gap equation (|48|) with p^l"^ ^ 1, we arrive at the 
following equation: 

ln^-~lnA2£2, (55) 



p2p ^p2p 
i.e. 

2 2 , ,ln[(lnA2£2)/vrl ^ ^ 

What is the physical meaning of this relation? Let us recall that at g = gc, the 
NJL model is equivalent to the (renormalizable) Yukawa model In leading 
order in 1/Nc, the renormalized Yukawa coupling ay "* = [gy '')^/(47r), corre- 
sponding to the scale /i = is cty ^ = ir/ In A^£^ 0. Therefore the expression 
(B^) for the dynamical mass can be rewritten as 



, a 



^dyn - |e5|— — In ^pxT • (57) 
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Thus, as has to be the case in a renormahzable theory, the cutoff A is removed 
through the renormahzation of parameters (the couphng constant, in this case) 
from the observable mdyn. 



Let us now consider the case of supercritical values of g {g > gc). In this case an 
analytic expression for mdyn can be obtained for a weak magnetic field, satisfying 
the condition \eB\^/'^ /m}-^^^ <^ 1, where rrS^^^ is the solution to the gap equation 



"dyr 

)l|) with B = Q. Then we find from Eq. (^Sf): 



2 / (0) \2 



\eB\ 



3(<l)Mn(A/m(5j2 



(58) 



i.e. mdyn increases with \B\. Q Notice that in the near-critical region, with 
9 — 1, this expression for m\^^ can be rewritten as 



2 ( (0) \2 

^dyn — l^^dynJ 



l + ^a'-^-^ lei? 



dyn; 



(59) 



where, to leading order in 1/A^c, Oy"'^^"^ = vr/ ln(A/mdyn)^ — ^/ lii(A/mdyn)^ is 
the renormalized Yukawa coupling related to the scale /i = mdyn- 

In conclusion, let us discuss the validity of the LLL dominance approximation in 
more detail. As Eq. ( ^3|) implies, the dynamical mass is 

mdyn = (0|a|0) = -G(0|V^^|0) . (60) 

Calculating the condensate (O|'i/''0|O) in the LLL dominance approximation, we 
find (see Eqs. (|I|) and (|19|)): 



,mdyn|eS| A^ 

Myn 

i.e. 



^"dyn = N^G — — In , (61) 
A-K^ mi 



m^y„ = A^exp(- , ) . (62) 



^ '\ \eB\N,G 

Comparing this relation with Eq. (|52|) , corresponding to dynamics with g outside 
the near-critical region, we conclude that this approximation reproduces correctly 
the essential singularity in mjyn at G = 0. On the other hand, the coefficient of 
the exponent in Eq. ([5^) gets a contribution from all the Landau levels. 



•^The fact that in the supercritical phase of the NJL model, md™ increases with \B\ has 
been already pointed out by several authors (for a review, see Ref. |l^). 
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In the near-critical region, with 1—g ^ In A^£^, all the Landau levels become 
equally important (see Eq. (^)), and the LLL dominance approximation ceases 
to be valid. This also happens (in a weak magnetic field) at supercritical values 
of g, g > gc, when the dynamical mass is generated even without a magnetic field. 

Thus we conclude that, like the BCS approximation in the theory of superconduc- 
tivity, the LLL dominance approximation is appropriate for the description of the 
dynamics of weakly interacting fermions. This point will be especially important 
in QED, considered in Section 9. 
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6. THE NJL MODEL IN A MAGNETIC FIELD: 
THE KINETIC TERM 



Let us now consider the kinetic term in the effective action (Bl). The chiral 
?7l(1) X f^R(l) symmetry imphes that the general form of the kinetic term is 



ipjd^Pj){pid^Pi) 



(63) 



where p = (cr, vr) and F^^ ^Fl^^ are functions of p^. We found these functions by 
using the method of Ref. [T^l. The derivation of Ck is considered in Appendix 
A. Here we shall present the final results. 



The functions F{ 



F^" are F[ 



^11 



p22 



n22 



gi^^F!^^, Fr = g^^'F^^ with 



87r2 



In ip 



.2«2 



+ 1 



1 1 

7^ ~ ^ ^ 3 



1 

p2 



1 



247r2 



P 



87r2 



.21)2 



C 2, 



o2/)2 



P 



2 02 



+ 1 



.2/)2 



47r 



- pH^ - pH^ + 1 



(64) 



where i\}[x) = d{\nr{x))/ dx. (We recall that the magnetic length is 



\eB\-^'\] 



As follows from Eqs. (^) and (|6^) , the propagators of a and vr in leading order 
in 1/Nc have a genuine (3 + l)-dimensional form. This agrees with the general 
arguments in support of the (3 + 1) -dimensional form of propagators of neutral 
particles in a magnetic field considered in Section 3. We shall see in Appendix 
B that this point is crucial for providing the reliability of the 1/Ac expansion in 
this problem. 

Now, knowing the effective potential and the kinetic term, we can find the dis- 
persion law for the collective excitations a and vr. 
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7. THE NJL MODEL IN A MAGNETIC FIELD: THE SPECTRUM 
OF THE COLLECTIVE EXCITATIONS 



The derivative expansion for the effective action is rehable in the infrared region, 
with k < Ft,, where ~ mdyn is the decay constant of vr. Therefore it is 
appropriate for the description of the low-energy dynamics of the gapless NG 
mode TT, but may, at most, aspire to describe only qualitatively the dynamics of 
the a mode. Nevertheless, for completeness, we will derive the dispersion law 
both for TT and a. 

We begin by considering the spectrum of the collective excitations in the subcrit- 
ical, g < Qc, region. At g <^ Qc = 1 "we find the following dispersion laws from 
Eqs. 



3D and 
E„ : 



m 



dyn 



eB 



In 



\eB\ 



12 m, 



dyn 



kl 



1/2 



\eB\ 



dyn 



1/2 



(65) 

(66) 



with mdyn defined in Eq. (|52[). Thus the vr is a gapless NG mode. Taking into 
account Eq. (0), we find that the transverse velocity |v_l| = \dET,/d]i.±\ of vr is 
less than 1. 

In the near-critical region, with g — 1 (where the NJL model is equivalent 
to the Yukawa model), we find from Eqs. (P7|), (|F7|), ( p5| ) and 



Et, 



E. 



1 



4"^dyn 



Inn /ay 
1 + 



ki +k 



3 Invr/ay 



(£-1) 



1/2 



+ 1- 



3 Invr/ay 



k^ +k:: 



1/2 



(67) 
(68) 



Thus, as has to be the case in a renormalizable model, the cutoff A disappears 
from the observables E^, and E^^. Note also that the transverse velocity v± of vr 
is again less than 1. 

Since the derivative expansion is valid only at small values of k±, k^, satisfying 
k±, k^ < mdyn, the relations (|65D-(|68D are valid only in the infrared region (as 
we already indicated above, the relations (|66D and ( |68D for a are only estimates). 
In particular, since, as i? ^ 0, the dynamical mass mdyn and the decay constant 
Fn ~ mdyn go to zero in the subcritical region, these relations cease to be valid, 
even in the infrared region, at B = 0. 

Let us now turn to the supercritical region, with g > gc. The analytic expressions 
for Ft, and E^^ can be obtained for weak magnetic fields, satisfying lei?!^/^ <^ 
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^dyn' where m^yn dynamical mass of fermions at S = 0. We find from Eqs. 

(0), (0) and ~ 



En 

E, 



1 - 



{eBf 



Am 



3m^y^ ln(A/mdyn' 
1 



k +k 



1/2 



dyn 



1 + 



1 



31n(A/mdyn 
ll{eBf 



2 1-^^ 

2 ^ 9m^y^ln(A/mdyn)2 



45m^yj^ ln(A/mdyn)^ 



nl/2 



(69) 



(70) 



with mdyn given in Eq. (|^). One can see that the magnetic field reduces the 
transverse velocity v± of rr in this well. 

Since in the supercritical region, with g — Qc 1, the renormalized Yukawa 

/ ln(A/mdyn)^ to leading order in Eqs. ( |69| ) and (^) 



coupling is a^'^^"'' = tt 
can be rewritten in that region as 



E-n 

E„ 



1 ^,(>»dyn) jeB) 

1 Cty 



4"^dyn 



("Idyii) 



k' +k 



1/2 



(71) 



1 + a'y '''''' ( — + 



A{eBf 



371 Ovrm^yn 



+ 1 - a 



(72) 



Let us now return to the discussion of the connection between the dynamics in 
the (3 + l)-dimensional NJL model in a magnetic field (with a weak coupling 
g <^ 1) and that in the Gross-Neveu (GN) model. Comparison of Eqs.(^) 
and (p3D suggests that the infrared dynamics in these two models may be sim- 
ilar. The GN model is asymptotically free, with the bare coupling constant 
(^(0) = 27i/A^c ln(A2/mdyn) — as A ^ oo. Therefore there is the dimensional 
transmutation in the model: in the scaling region, at -C 1, the infrared 
dynamics with (InA^/A;^)"^ -C 1 is essentially independent either of the cou- 
pling constant G^^^ or the cutoff A; the only relevant parameter is the dynamical 
mass mdyn (which is an analogue of the parameter Aqcd in QCD). One might 
expect that a similar dimensional transmutation should take place in the (3 + 1)- 
dimensional NJL model in a magnetic field: at \eB\GNc <C 1 (see Eq . (|52|) ) , the 
infrared dynamics should be essentially independent of the magnetic field and the 
coupling constant G. 

The situation in this model is however more subtle. As follows from Eq.([T8|) (with 
m replaced by mdyn), the fermion propagator in the infrared region is indeed 
independent (up to power corrections ~ (k5^/|ei?|)", n > 1) of the magnetic 
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field. However, the dependence of tlie propagator of the NG mode vr on \eB\ is 
essential. Though, as follows from Eqs.(p2l) and (|65D, the maximum transverse 



velocity |v_|_| = ^^J{m'^^^/\eB\)ln{\eB\/7rm'^y^) of tt is small at \eB\GNc -C 1, it is 
crucial that it is nonzero. The latter provides the (3 + 1) -dimensional character 
of the NG propagator and, as was already explained in Sec. 3, this in turn is 
crucial for the realization of spontaneous chiral symmetry breaking in the model. 

We recall that in the (l + l)-dimensional GN model with the chiral Ul{1) x Ur{1) 
symmetry, there is no spontaneous chiral symmetry breaking, despite the fact that 
the fermion dynamical mass is nonzero: in accordance with the MWC theorem 



T^ , there are no NG bosons in this model. As was suggested by Witten (see the 



second paper in Ref.|T^), in the GN model, the Berezinsky-Kosterlitz-Thouless 



(BKT) phase is presumably realized, with the NG bosons being replaced by the 
BKT vortex collective excitations. 

Therefore, in (3 + l)-dimensional NJL model, the magnetic length £ = \eB\~^^'^ 
acts as a (physical) regulator: at finite i, there is spontaneous chiral symmetry 
breaking, the propagator of NG mode vr has a (3 + l)-dimensional form and, as a 
result, the 1/Nc expansion is reliable (see Appendix B). In a sense, the magnetic 
length i plays here the same role as the e parameter in the (2 + e)-dimensional 
GN model [[gg. 



It is noticeable that the observables E^, and E^j in Eqs. (^5|) and (|6^) do not depend 
explicitly on cutoff A. It happens because, due to the relation (|52D for mdyn, the 
term l/p^£^|p2^p2=^2 dominates in the functions = Fp and ^2°° = Fg^^ 

in the kinetic term (see Eq. (|6^) ). In this dynamical regime, there is a strong 



hierachy of two scales: m'^y^ <C In a sense, with respect to the infrared 

dynamics with k"^ -C \eB\, the scale \eB\ plays the role of a (physical) ultraviolet 
cutoff. 

Let us consider the renormalization group properties of this dynamics. 

The reliability of the 1/Nc expansion implies that the infrared dynamics in this 
model is (at least at sufficiently large N^) under control: the effects connected 
with nonleading orders in l/N^ are small (see Appendix B). Let us consider 
the running Yukawa coupling in this model. In leading order in 1/Nc, the bare 
Yukawa coupling and both the renormalization constant of the fermion field and 
that of the Yukawa vertex equal one (see Appendix B). Therefore the behavior 
of the running Yukawa coupling is determined by the renormalization constant 

^Comparing the effective actions in the (3 + l)-dimentsionai NJL modci in a magnetic field 
and (1 + l)~dimentsional GN model, one can show (at least formally) that, as \eB\ oo, the 
NJL model, with Nc colors becomes equivalent to the GN model with the number of colors 
Nc = NcN, where N = \eB\Si2/2TT ^ cjo is the number of states at the LLL {S12 is the square 
in the a;ia;2-plane). This interesting limit will be considered in more detail elsewhere. 
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of the chiral field p: ~ Eq.® implies that Z^'^) ' ~ A/'e|e5|/87rV^ 

at /i ~ JTidyn- Therefore the running Yukawa coupling is weak in the infrared 

region: gfy "* ~ ^Js^l'^^'^ /Nc\eB\ and gy"^ — > as |e-B| (and therefore also A) goes 
to infinity, i.e. there is a Gaussian infrared fixed point. 

Examining the effective potential and the kinetic term in the effective action, 
one can show that a similar behavior occurs for the rest of running couplings, 
describing self-interactions of the chiral field p. 

We emphasize that this discussion relates only to the NJL model with a weak 
coupling constant. In the case of the NJL model with a near-critical g, the 
situation is different. In that case, one finds from (p^ ) that the renormalization 
constant Z^'^^ is Z^J^^ ~ (A''c/47r^) ln(A/yu). Therefore the behavior of the 
running coupling is similar to that in the NJL model without magnetic field (and 
with the same near-critical g) P]: g'^^ ~ 2n / ^N^AnX/Jl, i.e. there is the usual 
Gaussian infrared fixed point in this case. 

As was already indicated in the previous section, the difference between these 
two dynamical regimes is connected with the point that, while a.t g <^ 1 the LLL 
dominates, at a near-critical g, all Landau levels become relevant. 

In the next section, we will derive the Bethe-Salpeter equation for the NG mode 
TT in the LLL dominance approximation. This will yield a further insight into 
the physics of the dimensional reduction D D — 2 in the dynamics of fermion 
pairing in a magnetic field. 
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8. THE NJL MODEL IN A MAGNETIC FIELD: 
THE BETHE-SALPETER EQUATION 
FOR THE NG MODE 



The homogeneous Bethe-Salpeter (BS) equation for the NG bound state vr takes 
the form (for a review, see Ref. P]): 

XAB{x,y;P) = -i j d'^xid^yid'^X2d^y2GAAAx,xi)KA^Bv,A2B2{xiyi,X2y2) 

X XA2B2{.X2,y2]P)GB^B{,yuy) , (73) 

where the BS wave function XAB{x,y; P) = {0\TipA{x)iljB{y)\P',T!') and the fer- 
mion propagator GAB{x,y) = (0|T'?/^^(x)'?/^b(?/)|0); the indices A = (na) and 
B = {m(3) include both Dirac (n, m) and color (a, (3) indices. Notice that though 
the external field A^^^ breaks the conventional translation invariance, the total 
momentum P is a good, conserved, quantum number for neutral bound states, in 
particular for the vr (see Section 3). Henceforth we will use the symmetric gauge 
(0). 

In leading order in l/Nc-, the BS kernel in the NJL model is 0: 

KAiBi;A2B2{Xiyi,X2y2) = G[5ajBi5b2A2 + 5ai/3i5/32a2(n5)nimi(«75)m2n2 

— SaiA2^B2Bi — <^aia2'^/32/3i(^75)nin2(n5)m2mi] 

X 6\x^ - yi)5\x^ - X2)5\x^ - t/s) • (74) 

Also, as was already indicated in Section 5, in this approximation, the full-fermion 
propagator coincides with the propagator S (|^) of a free fermion (with m = 
mdyn) in a magnetic field. 



Let us now factorize (as in Eq. ([T3|) ) the Schwinger phase factor in the BS wave 
function: 

Xab{R, r; P) = exp(^er'^A^^*(i?))e"^'^^Xnn^(^, r; P) , (75) 

where we introduced the relative coordinate, r = x — y, and the center of mass 
coordinate, R = {x + y)/2. Then Eq. ([7^) can be rewritten as 



XnM r; P) = -iN,G J d^R.Snm {^- + R-Ri 



■'nimi 



tr(x(Pi,0;P)) 



(75)nimitr(75x(Pi,0;P)) - ^Xnimi(i?l,0;P) 



+ Jj-{l5)nrn2Xn2m2{Rl,0;P){l5)m2m^ 



c 

x^^i™(^ - P + Pi) exp[-^erM^^^*(^ - ^i)] exp[iP(P - P^] . (76) 
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The important fact is that the effect of translation symmetry breaking by the 
magnetic field is factorized in the Schwinger phase factor in Eq. (fTSD, and Eq. 
(^) admits a translation invariant solution: Xnm{R,r; P) = Xnm{r,P). ^ Then, 
transforming this equation into momentum space, we get: 



d^q±d'^R±d k±d k^l 



X exp[i(Px - q^)R^]5„„, + ^, p± + eA*^"*(Rx) + ^ 

^nimi tr(x(A;;P)) - (75)nimi tr(75x(x; P)) - ■^Xmmi{k;P) 



X ^^,™(p||-5l,p^ + eA^^*(Rx)-^^ (77) 



where = {p ,p ), p_|_ = (p ,P )■ Henceforth we will consider the equation with 
the total momentum — > 0. 

We shall consider the case of weakly interacting fermions, when the LLL pole ap- 
proximation for the fermion propagator is justified. Henceforth, for concreteness, 
we will consider the case eB > 0. Then 

S{p) ^^eM-<^'pl)^^^il - ^iV) (78) 

P\\ " ''^dyn 

(see Eq. (p!8|)), where = p'^^O— p^^s^ ^^^^ transforms into the following 

one: 



= ^e-^^Pi / d^A^d-k^d\,e-''^\x - .7S^) 
X P[p(A;||,k^)](l-z7S') , (79) 



where 

P(P|| , P±) = (P|| - "^dyn)x(P|| , P±) (P|| - "^dyn) (80) 

with x(P||,p±) = x(P||) P±j -P)lp=o; operator symbol P[p] means: 

mh-^.)\ = t<|l^p(*ii.kx)||i^) 



^This point is intimately connected with the fact that these bound states are neutral: the 
Schwinger phase factor is universal for neutral fermion-antifermion bound states. 
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, k\\ +mdyn , . +mdyn\ 



1 +mdyn , +mdy„ 



1 fcy+mdyn , ^A;||+mdyn 
+ ]v75 / p(fc||,k^) " / 75- (81) 

Eq. (|79D implies that p{p\\,p±) = exp(— £^p5_)v?(p||), where (p{p\\) satisfies the 
equation: 

¥^(^11) = 1^ / rf^A;||(l -.7V)^b(A:||)](l - ) • (82) 

Thus the BS equation has been reduced to a two-dimensional integral equation. 
Of course, this fact refiects the two-dimensional character of the dynamics of the 
LLL, that can be explicitly read off from Eq. (|7Sp. 

Henceforth we will use Euclidean space with ^4 = — In order to define the 
matrix structure of the wavefunction ip{p\\) of the vr, note that in a magnetic field, 
in the symmetric gauge (10), there is the symmetry 5*0(2) x 5*0(2) x V, where 



the 5*0(2) X 5*0(2) is connected with rotations in the xi — X2 and X3 — X4 planes 
and V is the inversion transformation X3 —X3 (under which a fermion field 
transforms as ip ^ ^7573"^) • This symmetry implies that the function if{p\\) takes 
the form: 

'f{P\\) = 75(^ + ^7172^ + P\\C + i7i72]5||D) (83) 

where p\\ = P373+P474 and A, B, C and D are functions of (7^ are antihermitian 
in Euclidean space). Substituting expansion (^) into Eq. (^), we find that 
B = —A, C = D = 0, i.e. f{p\\) = ^75(1 — ■^7172)- Q The function A satisfies the 
equation 

The solution to this equation is A{p) = constant, and introducing the ultraviolet 
cutoff A, we get the gap equation for m 



2 . 

dyn- 



N^G r^' dk 



A'k'^P Jo fc2 + m^y„ 



^5) 



^The occurrence of the projection operator (1 — ij^j'^)/2 in the wave function ip reflects the 
fact that the spin of fermions at the LLL is polarized along the magnetic field. 
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It leads to the expression 



for m\y^ which coincides with expression ( p^D derived in Section 5. 

The integral equation (|8^ ) can be rewritten in the form of a two-dimensional 
Schrodinger equation with a 5-like potential. Indeed, introducing the wave func- 
tion 



we find 



where 



N C 

-A + ml^--^6l{r))^ir)=0 (87) 



and 



A (27r) 

Notice that in the same way, one can show that in the (2 + l)-dimensional NJL 
model , the analog of Eq. (^) has the form of a one-dimensional Schrodinger 
equation: 

(-^ + <n - -^^K[r)Y{r) = (89) 



with (5a (r) = J\ l^e It leads to the following gap equation for rri^y^. 

NrG dk 



(90) 



which yields expression (^) for m^y^ i^i the limit A — >• cxo: m^^^ = (|ei?|A''cG'/2 

Thus, since (— "^dyn) pl^-ys the role of the energy E in the Schrodinger equation 
with a negative, i.e. attractive, potential, the problem has been reduced to finding 
the spectrum of bound states (with E = — m^y^ < 0) of the two-dimensional 
and one-dimensional Schrodinger equation with a short-range, (5-like, potential. 
(Actually, since only the largest value of m^yn defines the stable vacuum one 
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has to find the value of the energy for the lowest stationary level.) This allows us 
to use some general results proved in the literature that will in turn yield some 
additional insight into the dynamics of chiral symmetry breaking by a magnetic 
field. 

First, there is at least one bound state for the one-dimensional {d = 1) and 
two-dimensional {d = 2) Schrodinger equation with an attractive potential [0. 
Second, while the energy of the lowest level E{G) is an analytic function of the 
coupling constant, around G = 0, at ci = 1, it is non-analytic at G = at = 2 
0. Moreover, at d = 2 for short-range potentials, the energy E{G) = —m'^y^{G) 
takes the form E{G) ~ — exp[l / (aG)] (with a being a positive constant) as G ^ 

0- 

Thus the results obtained in the NJL model at D = 2 + 1 and D = 3 -\- 1 agree 
with these general results. 

The fact that the effect of spontaneous chiral symmetry breaking by a magnetic 
field is based on the dimensional reduction D — > D — 2 in a magnetic field 
suggests that this effect is general, and not restricted to the NJL model. In the 
next section, we shall consider this effect in (3 + 1) -dimensional QED. 
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9. SPONTANEOUS CHIRAL SYMMETRY BREAKING 
BY A MAGNETIC FIELD IN QED 



The dynamics of fermions in a constant magnetic field in QED was first considered 
by Schwinger In that classical work, while the interaction with the external 
magnetic field was considered in all orders in the coupling cosntant, the quantum 
dynamics was treated perturbatively. There is no dynamical chiral symmetry 
breaking in this approximation In this section we reconsider this problem, 
treating the QED dynamics non-perturbatively, and show that, in ladder and 
improved ladder approximations, a constant magnetic field leads to spontaneous 
chiral symmetry breaking. 

We will use the same strategy as in the previous section and derive the BS equa- 
tion for the NG modes. 

The Lagrangian density of massless QED in a magnetic field is: 

C = AF^'^F,^ + ^[^P,{^^^DM , (91) 
where the covariant derivative is 

D^ = d,- ze(A^^* + A^) , A^;' = (o, - |x2, o) , (92) 

i.e. we use the symmetric gauge (IC). Besides the Dirac index (n), the fermion 
field carries an additional flavor index a = 1,2,..., Nf. Then the Lagrangian 
density ( pl| ) is invariant under the chiral SUj,{Nf) x SUji{Nf) x ?7l_|_r(1) (we 
will not discuss the dynamics related to the anomalous, singlet, axial current Js^ 
in this paper). Since we consider the weak coupling phase of QED, there is no 
spontaneous chiral symmetry breaking at i? = [1^. We will show that the 
magnetic field changes the situation dramatically: at B ^ the chiral symmetry 
SULiNf) X SUniNf) breaks down to SUyiNf) = SUL+niNf). As a result, the 
dynamical mass mdyn is generated, and Nj — 1 gapless bosons, composed of 
fermions and antifermions, appear. 

The homogeneous BS equation for the Nj — 1 NG bound states takes the form 
i: 

Xab^^^V'^P) = ~^ / d'^xid'^yid'^X2d'^y2GAAiix,xi)KA^Bi;A2B2ixiyi,X2y2) 

X XMB2i.^'i,y2]P)GB^B{y2,y) , (93) 



where the BS wave function x^ab = {^\Ti^A{x)'il)B{y)\P'-, /3 = 1, . . . , Nj — 1, and 
the fermion propagator GAB{x,y) = {0\TipA{x)ipB{y)\0); the indices A = (na) 
and B = {mb) include both Dirac (n, m) and flavor (a, h) indices. 
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As will be shown below, the NG bosons are formed in the infrared region, where 
the QED coupling is weak. This seems suggest that the BS kernel in leading order 
in a is a reliable approximation. However, because of the (1 + 1) -dimensional 
form of the fermion propagator in the infrared region, there may be also relevant 
higher order contributions. We shall return to this problem below, but first, we 
shall analyze the BS equation with the kernel in leading order in a. 

The BS kernel in leading order in a is |0: 

X6{XI - X2)6{yi - y2) + ^^^ia^a^bA2a2ln^m^lm2n2^^^^>^i.^l " ^2) 

x5(xi - yi)5{x2-y2) , (94) 
where the photon propagator 

(A is a gauge parameter). The first term on the right-hand side of Eq. (0) 
corresponds to the ladder approximation. The second (annihilation) term does 
not contribute to the BS equation for NG bosons (this follows from the fact that, 
due to the Ward identities for axial currents, the BS equation for NG bosons can 
be reduced to the Schwinger-Dyson equation for the fermion propagator where 
there is no contribution of the annihilation term [^]). For this reason we shall 
omit this term in the following. Then the BS equation takes the form: 

P) = -4.a j d^M^^^i., ■^.K..,l^„xUM. v.: P) 

X ^b2bilrn2mi 

SBiB{yi,y)T>^„y{yi - Xi) , (96) 

where, since the lowest in a (ladder) approximation is used, the full fermion 
propagator Gab{x, y) is replaced by the propagator 5* of a free fermion (with the 
mass m = m^y^) in a magnetic field (see Eqs. (|13|), (|1^). 

Using the new variables, the center of mass coordinate R = ^-j^, and the relative 
coordinate r = x — y, Eq. (^) can be rewritten as 

XnmiR, r; P) = -47r« J d'Rid'nSnn, (^R-Ri + '-Y^) 

~ ( T — r\ \ 

(^l,'"i;^)7m2mi Sm^my^ — ^ -R + R\^D y^y[—r\) 

X exp[-ze(r + r^YA^^\R - R^)\ x exp[iP(i? - R^)\ . (97) 

Here S is defined in Eqs. (|13|) and (|ll), and the function Xnm{R-,'r \ P) is defined 
from the equation 

X^ABi^,y;P)^{0\TMx)^PB{y)\P,P) 

= Xie-'''^exp[zer^A^;\R)]xnmiR,r;P) (98) 
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where are Nj - 1 flavor matrices (tr(A^A'^) = 2S/3^; /?,7 = l,...,Nj -1). 
The important fact is that, like in the case of the BS equation for the vr in the 
NJL model considered in the previous section, the effect of translation symmetry 
breaking by the magnetic field is factorized in the Schwinger phase factor in 
Eq. (|98D, and Eq. (^) admits a translation invariant solution, Xnm{R,r; P) = 
x(r; P). Then, transforming this equation into momentum space, we get 



Xnm{p; P) = -47ra J 



(fq±(fR±(fk±(fkii 



X exp[^(Px - q±)R±]Snn, (p\\ + y . + eA'=^*(R^) + ^ 

/ P\\ Q-L 

X P^,(A;|| -p||,kx-p±-2eA^^*(R^)) (99) 

(we recall that p\\ = {p'^,p^), p± = (p^,P^))- Henceforth we will consider the 
equation with the total momentum 0. 



The crucial point for further analysis will be the assumption that mdyn << y |e-B| 
and that the region mostly responsible for generating the mass is the infrared 



region with k < m^y^ <^ y |ei?|. As we shall see, this assumption is self-consistent 

(see Eq. @). The assumption allows us to replace the propagator Snm in Eq. 
(P^) by the pole contribution of the LLL (see Eq. (^)), and Eq. (^9]) transforms 
into the following one: 



P{P\\,P±) = ^,e-''^-Jd'A^d'k^d%e-''^'i{l-^JV)Y 

X (1 - ^7V)V,Ah - Plh - A,) , (100) 

where p{p\\,p±) = (P|| -"^dyn)x(p)(P|| -"^dyn)- Eq. ( [lOOD implies that p(p||, p±) = 
exp(— £^p5_)v5(p||), where f{p\\) satisfies the equation 

X ^4^^nl-^lV)Dl{h-p,). (101) 



Here 



-Pll) = / d''k^exp(^-^^^V,Uk\\ -P\\,^±) ■ (102) 
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Thus, as in the NJL model in the previous section, the BS equation has been 
reduced to a two-dimensional integral equation. 

Henceforth we will use Euclidean space with /C4 = —ik^. Then, because of the 
symmetry S0{2) x SO (2) xP in a magnetic field, we arrive at the matrix structure 
for ipipii): 

'^{P\\) = 75(^4 + 271725 + p\\C + i'ji-f2P\\D) (103) 
where A,B,C and D are functions of . 

We begin the analysis of Eq. (|101|) by choosing the Feynman gauge (the general 
covariant gauge will be considered below). Then 

r°o dxexp{-fx/2) 



and, substituting the expression (|103|) for (p{p\\) into Eq. (|101|) , we find that 
B = -A, C = D = 0, i.e. 

(^(Pll) =^75(1- 27172) , (105) 
and the function A satisfies the equation: 

a f d?kA{k) dxexp{-fx/2) 



"^^'^ = 2^^ J k^T^Jo (k-pr + x ^'''^ 

(henceforth we will omit the symbol ||; for the explanation of physical reasons of 
the appearance of the projection operator (1 — 27i72)/2 in (p{p), see footnote P). 

Now introducing the function 

J {2tt)^ k^ + miy^ 
we get the two-dimensional Schrodinger equation from Eq. (|106|) 

{-A + ml^ + V{r))^{r)=0 , (108) 
where the potential V^(r) is 

T// \ " f j2 inr dxexp{-x/2) a , _^i2t^ (r r- 

^ 2 s ^ 2 



32 



{Kq is the Bessel function). The essential difference of the potential ( |109| ) with 
respect to the S-\ike potential ( pBD in the NJL model is that it is long range. 
Indeed, using the asymptotic relations for Ei{x) 0, we get: 

2a 1 



Vir) 
Vir) 



oo 



a / 



7 + In ■ 



2e 







;iio) 



where 7 ^ 0.577 is the Euler constant. Therefore, the theorem of Ref.[0 (asserting 
that, for short-range potentials, —E{a) = m'^ (a) ~ exp(— 1/aa), with a > 0, 



as a — i> 0) cannot be applied to this case. In order to find m'^y^{a), we shall use 
the integral equation ( |106| ). This equation is analyzed in Appendix C by using 
both analytical and numerical methods. The result is: 

TT / TT \ 



"^dyn = C \/\eB\ exp 



111 



2 \2aJ 

where the constant C = 0(1). Note that this result agrees with the analysis of 
Ref. |]19[ where the analytic properties of E{a) were studied for the Schrodinger 
equation with potentials having the asymptotics V^(r) const /r^ as r — 00. 



Since 



lim 



exp[—l / a^/a\ 



00 



;ii2) 



+0 exp[— l/a'a] 

at a, a' > 0, we see that the long-range character of the potential leads to the 
essential enhancement of the dynamical mass. 



Let us now turn to considering the general covariant gauge (|95f ). As is known, the 
ladder approximation is not gauge invariant. However, let us show that because 
the present effect is due to the infrared dynamics in QED, where the coupling 



constant is small, the leading term in ln(mjyj^£^), 
same in all covariant gauges. 



ln(m^ ^'^'^ 



-TT 



7r/2a;, is the 



Acting in the same way as before, we find that the wave function (p{p) now takes 
the form 



V.(p)=75(l-nV)(^(p')+pO(p2)) 
where the functions A{p'^) and C{p'^) satisfy the equations: 



(fkA{k^ 



a 



277^ J + m^yjj JO 



dx{l - Xx^/A) exp(-x£V2) 



Att'^J k^ + 



m 



dyn 



2e - (kp) 



(k — p)^ + X 
P(kp)- 



(113) 
(114) 



X 



dx exp(— x£^/2) 

[{k-py + xf 



(115) 
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One can see that the dominant contribution on the right-hand side of Eq. ( |1 14|) 
(proportional to [Inm^yj^f^jz g^^^ formed at small k'^) is independent of the gauge 

parameter A. Thus the leading contribution in ln(m^yn£^), \n{m\j^i'^) ~ —n^7r/2a, 
is indeed gauge-invariant (for more details see Appendix C). 



This concludes the derivation of Eqs. (|), ( |108[ ) and ( p.09| ) describing spontaneous 
chiral symmetry breaking by a magnetic field in ladder QED. 

Note the following important point. Because of the exponent exp{—i'^x/2) in 
integral equations ( p.06| ) , ( |114|) and (|115|) , the present effect is connected with the 
infrared dynamics in QED: the natural cutoff in this problem is \eB\. The fact 
that the non-pertubative infrared dynamics in this problem decouples from the 
ultraviolet dynamics is reflected also in the asymptotic behavior of the functions 
y4(p2) and C(p^): as follows from Eqs. (|106|) , ( 114 ) and (|115|) , these functions 



rapidly {A{p ) ~ 1/p , C{p ) ~ A/p ) decrease as p oo. Therefore, since 
at p^ ^ |ei?| the magnetic field essentially does not affect the behavior of the 
running coupling in QED the coupling constant a in equations ( |106D , ( |114D 
and (|115D has to be interpreted as the value of the running coupling related to 
the scale fi^ ~ \eB\. 

Up to now we have considered the ladder approximation in QED in a magnetic 
field. As was already mentioned above, because of the (l + l)-dimensional form of 
the fermion propagator in the infrared region, there may be also relevant higher 
order contributions. In particular, the dimensional reduction may essentially af- 
fect the vacuum polarization, thus changing the behavior of the running coupling 
in the infrared region. We will show that this is indeed the case. Actually, in 
this problem, because a magnetic field breaks Lorentz and 5*0(3) rotational sym- 
metries, it is more appropriate to speak not about a single running coupling but 



about a running coupling tensor, 
the magnetic field. 



i.e. about the full photon propagator Vn^, in 



Let us consider the photon propagator in a strong magnetic field in the infrared 
region {\eB\ ^ w^dyn; l^yl; l^ll)' with the polarization operator calculated in 
one-loop approximation 



21 1 . One can rewrite it in the following form: 



V 



U{kl 



+ 



1 



1 



p + Hn(H) ^^^^^ 



a 



\eB\ 



m2y^ 



fc2 



k\\kl 



A k^ki, 
¥~k^^ 

4 



6m 



dyn 



4"^dyn^y 



X 



;ii6) 
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X In 



{kfy - Aml^^q + kj 



;ii7) 



where the symbols ± and || in (yf^j, are related to the (1, 2) and (0, 3) components, 
respectively. The asymptotic behavior of Il{kf) is: 



ml) 
ml) 



a \eB\ 
37rm^y^ 
2a \eB\ 



at |fc|| I <^ ^dyn? 



at 1^1 1 > "^dyn- 



;ii8) 
;ii9) 



Note the following characteristic points: 



correspond to the one-loop contribution with 
117|) is the leading term in the 



a) The expressions ( |116|) , 

the fermions from the LLL. Actually, expression 
l/|ei?| expansion of the one-loop polarization operator. 

b) The polarization effects are absent in the transverse components of 1^^,^. This 
is because the spin of fermions from the LLL is polarized along the magnetic field 
(see footnote ^). Indeed, this implies that in the QED-vertex, with two fermions 
from the LLL, the photon spin equals zero along the magnetic field, i.e. only the 
longitudinal, (0,3), components, are relevant in this case. 

c) There is a strong screening effect in the ((?l[^ — fc|l/c|/A;y)-component of the 
photon propagator. In particular there is a pole at kj = in n(fc|) as m^^^ — > 0: 
this is a reminiscence of the Higgs effect in the (1 + l)-dimensionai massless QED 
(Schwinger model) 12^ 



d) Is the dynamics in QED in a magnetic field similar to that in the Schwinger 
model, as \eB\ —>■ oo7 The answer to this question is "no". The crucial point 
is that, unlike the NJL model, there is a neutral fundamental field in the QED 
Lagrangian: the photon field. As one can see from Eq.( p.l6| ), the photon propa- 
gator does not have a (1 + l)-dimensional form even as \eB\ — >■ oo: there is the 



quantity k 



k"^ 



kj_ — k^ in the denominator of the photon propagator. This 



leads to the interaction which is very different from that in the Schwinger model. 
Notice that the interaction at long distances in this problem is much weaker than 
that in the (1 + 1) -dimensional Schwinger model: while here, in Euclidean space, 
the potential V{r) is V{r) ~ — as r ^ cxd (see Eg. ( 110 )), it is V{r) ~ Inr in 
the Schwinger model. Thus the infrared dynamics in QED in a magnetic field is 
a complex mixture of (3 + 1) and (1 + 1) -dimensional dynamics. 

Let us now consider the BS equation in the so called improved ladder approx- 
imation: in this approximation, the free photon propagator is replaced by the 



one-loop propagator (|116|). As was already indicated above, the transverse com- 
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ponent of the photon propagator decouple from the fermions at the LLL. The 
longitudinal components lead to the following equation (in the Feynman gauge): 



a 



J ki + mi.„ J 



1 



X + 



+ 



+ 



X + 



+ (fc||-p||)2n£;(A;||-. 



(120) 



where II^; is the polarization operator ( |117|) in Euclidean space (compare with 
Eq. (|106D ) . The first term in the square brackets on the right-hand side of Eq. ( |120| ) 
comes from the (unscreened) /c|[/c||//c^-component; the second term comes from 
the (screened) {g^^^, — A;|[fcjj/fcy)-component. It is not difficult to show that at 
= the latter does not contribute to the dominant term (proportional to 
[Inm^yn^^]^) on the right-hand side of this equation. Indeed, one finds 



AiO) 



a 

An' 



^ + ^dyn 



dxe-'"^^^ 



+ 



x + kj X + kd + k?,UE{k\\) 



1211 



and, matching asymptotics ( |118|) , ( 119 ) at k 



the following estimate for the second term: 

d^h 



6m^yn in Euclidean space, we get 



a 
4^ 



'J kl + m%J 



6mj 

dyn 



— A(0) j dxe-'"^/' J 



dy 



X + kf^ + kjUEikw) 
dy 



(y + ^dyn)[^ + 2/(1 + a/STTmLn^^)] 



+ 



+ 



6m 



iy + ml^J{x + y + 2a/ne) 



dyn 



7"^dyn^ 



In^. 

a 



(122) 



Therefore the dominant contribution in this equation is connected with the un- 
screened fc|[/c|/A;^fc^ component (since is different from fc^, this term is not a 
gauge artifact wich can be removed by a gauge transformation). As a result, we 
get the same equation as in the ladder approximation but with a replaced by 
a/2. The expression for mayn is now given by Eq.( |lll| ) with a — > a/2. It is 
gauge invariant. Also, since a is the value of the running coupling at a definite 
scale, /i^ ~ this expression is renormalization group invariant. 

The present consideration shows that, despite the smallness of a, the expansion 
in a is broken in the infrared region in this model. There are two reasons for 
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that: the (1 + 1) -dimensional character of the fermion propagator in the infrared 
region and the smallness of the dynamical mass mdyn as compared to |ei?|^/^ 
{H\eB\/m%^) ~ 

It is well known that in massive QED, infrared singularities (on mass shell) in 
Green's functions can be completely factorized if a certain set of diagrams is 
summed up p|. The situation in massless QED is much more complicated: the 
set of relevant diagrams becomes essentially larger and, up to now, the problem 
of the complete description of the infrared singularities in that model remains 
unsolved. The smallness of mdyn in QED in a magnetic field makes this model 
closer to massless QED than to massive one. This implies that, despite the 
smallness of the coupling constant a, the infrared dynamics in this model is quite 
complicated. 

It is a challenge to define the class of all those diagrams in QED in a magnetic 
field which give a relevant contribution in this problem. Since the QED coupling 
constant is weak in the infrared region, this problem, though hard, seems not to 
be hopeless. 
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10. CONCLUSION 



In this paper we showed that a constant magnetic field is a strong catalyst of 
spontaneous chiral symmetry breaking in 3+1 and 2+1 dimensions, leading to the 
generation of a fermion dynamical mass even at the weakest attractive interaction 
between fermions. The essence of this effect is the dimensional reduction D 
D — 2 in the dynamics of fermion pairing in a magnetic field. In particular, the 
dynamics of NG modes, connected with spontaneous chiral symmetry breaking by 
a magnetic field is described by the two (one)-dimensional Schrodinger equation 
at D = 3 + 1 (D = 2 + 1): 

i-A + ml, + V{r))<iJ{r)=0, (123) 

where the attractive potential is model dependent and it defines the form of the 
dynamical mass as a function of the coupling constant. Since the general theorem 
of Ref. assures the existence of at least one bound state for the two- and one- 
dimensional Schrodinger equations with an attractive potential, the generation of 
the dynamical mass m^y^ takes place even at the weakest attractive interaction 
between fermions at D = 3 + 1 and D = 2 + 1. This general effect was illustrated 
in the NJL model and QED. 

In this paper, we considered the dynamics in the presence of a constant magnetic 
field only. It would be interesting to extend this analysis to the case of inhomo- 
geneous electromagnetic fields. In connection with this, note that in 2+1 dimen- 
sions, the present effect is intimately connected with the fact that the massless 
Dirac equation in a constant magnetic field admits an infinite number of normal- 
ized solutions with E = (zero modes) 0. More precisely, the density of the 
zero modes 



z>o = lim S-^N{E) 



s 



E=0 



(where is a two-dimensional volume of the system) is finite. As has been 
already pointed out (see the second paper in Ref. and Ref. |^), spontaneous 



flavor (chiral) symmetry breaking in 2+1 dimensions should be catalysed by all 
stationary (i.e. independent of time) field configurations with being finite. On 
the other hand, as we saw in this paper, the density 



i>Q = lim V 



.idN{E) 



V^oo dE 



E=0 



of the states with E = (from a continuous spectrum) plays the crucial role in 
the catalysis of chiral symmetry breaking in 3+1 dimensions. One may expect 
that the density z/q should play an important role also in the case of (stationary) 
inhomogeneous configurations in 3+1 dimensions. 
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As a first step in studying this problem, it would be important to extend the 
Schwinger results to inhomogeneous field configurations. Interesting results in 



this direction have been recently obtained in Ref. ||24 



In conclusion, let us discuss possible applications of this effect. 

Since (2 + 1) -dimensional relativistic field theories may serve as effective theories 
for the description of long wavelength excitations in planar condensed matter 
systems ||2^, this effect may be relevant for such systems. It would be also inter- 
esting to take into account this effect in studying the possibility of the generation 
of a magnetic field in the vacuum, i.e. spontaneous breakdown of the Lorentz 
symmetry, in (2 + 1) -dimensional QED |]26[ . 



In 3+1 dimensions, one potential application of the effect can be connected with 
the possibility of the existence of very strong magnetic fields {B ~ lO^^G) during 



the electroweak phase transition in the early universe ||2^. As the results obtained 
in this paper suggest, such fields might essentially change the character of the 
electroweak phase transition. 

Another application of this effect can be connected with the role of chromo- 
magnetic backgrounds as models for the QCD vacuum (the Copenhagen vacuum 



Yet another potentially interesting application is the interpretation of the results 
of the GSI heavy-ion scattering experiments in which narrow peaks are seen in 
the energy spectra of emitted e~^e~ pairs |2^. One proposed explanation |Q is 



that a strong electromagnetic field, created by the heavy ions, induces a phase 
transition in QED to a phase with spontaneous chiral symmetry breaking and the 
observed peaks are due to the decay of positronium-like states in the phase. The 
catalysis of chiral symmetry breaking by a magnetic field in QED, studied in this 
paper, can serve as a toy example of such a phenomenon. In order to get a more 
realistic model, it would be interesting to extend this analysis to non-constant 



background fields |31 



We believe that the effect of the dimensional reduction by external field configura- 
tions may be quite general and relevant for different non-perturbative phenomena. 
It deserves further study. 
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APPENDIX A 



In this Appendix, we derive the kinetic term Ck in the effective action F 
the NJL model. 



m 



The structure of the kinetic term 



is: 



(124) 



where p = (a, vr) and Fj"^, Fl^^ depend on the Ui,{l) x ?7R(l)-invariant p^ = ct^+tt^. 
The definition F = / d^xC and Eq. ( |124| ) imply that the functions -Fi"', -^2''^ 
determined from the equations: 



6cr{x)6a{0) 



6n{x)6Tr{0) 



a=const 
7r=0 



cr=const 
7r=0 



a=const 

TC=0 



a=const 
7r=0 



■d^dy5'^{x) 



(125) 
(126) 



Here F^ is the part of the effective action containing terms with two derivatives. 
Eq. (|3^) implies that F^ = F^. Therefore we find from Eq. (|126D that 



fii/ 



d'^xx'^x'^ 



57r(x)(57r(0) 



d'^xx^x" 



5^T 



57r(x)57r(0) 



(127) 



(henceforth we shall not write explicitly the condition cr=const., vr = 0). Taking 
into account the definition of the fermion propagator, 



iS-^ = iD-a 



:i28) 



we find from Eq. 



6^T 



6tt{x)6tt{0) 



-I tr(5(x, 0)i7^S(0,x)i7^ 

-i ii{S{x)i^^S{-x)i^^) 
d'^kd'^q 



i2ny 



e"^''tT{S{k)iYS{k + q)tY) 



(129) 



(the functions S{x) and S{k) are given in Eqs. (|T3]) and (|T^) with m = a). 
Therefore 



d^k 



(130) 
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In the same way, we find that 



Mm^-m^,^^^,^]. (131) 



' A J (27r)4 V ' 'dk^dK dk^dK 

Taking into account the expression for S{k) in Eq. (|1^) (with m = a), we get: 

d''S{k) 
dkodko 



6^S{k) 
dk^dks 



d^Sjk) 
dkjdkj 



2iN,t / dt ■ te^(*) [(2^^(^A;°)'[A;°7° - ^-f^ + a] 
Jo 

+3fcV - fcV + (^)(1 + Vl^l^T) - kx7x(l + 2i{ik'^)h) 

x{l + T')], (132) 

-2iN,t / dt ■ [(-2^^ (£A;')' [fc°7° - + a] 
Jo 

+fc°7° - 3A;V + f^)(l + Vl^l^T) - kx7x(l - 2i(£A:3)2t) 

x(1 + T2)], (133) 

-2zAr^£^ / dt ■ Te-^W[(-2iT(£A;-')2[fcV - '^V + 
Jo 

+fc0^o _ ^3^3 ^ ^ r/7S'T) - k^7^(l - 2i{ik^fT){l + T^) 
-2fcV(l + T2)], (134) 



where T = tant, r] = sign(e5), R{t) = UH{{k^f - (P)^ - a"^) - ifk^T, and 
j = 1,2 (there is no summation over j). 

Eqs. (|1|), ( [130D , and ([ml) imply that non-diagonal terms F/"^ and Fg"^ are equal 
to zero. Below, we will consider in detail the calculation of the function 
other functions Ff^, ^2^*" can similarly be found. 

Substituting expressions (0) and (|132|) into Eq. ( |130|) for we get 

^ ^ I d'k r r dTdte^^^'^+^^^'>hA[2z{ik')H{{ky - {k'f - a^) 



A'k'^ j Jo Jo 

+ 3{k^f - (P)2 - a^) - Bklil + 2i{ik'>)H)] , (135) 

where A = 1- TT, B = {1 + T'^){1 + T^), T = tanr. After integrating over k, 
we find: 

^00 _ roo roo drdt ^^^_,,2,2(,^,) 



X 



47r2yo Jo (r + t)2 



(136) 



Changing t and r to new variables. 



t='-{l + u), T='-{l-u), (137) 
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and then introducing the imaginary (dimensionless) proper time, s —is, we 
arrive at the expression: 

+ 





ds 


247r2 


Jo s 




r°° ds 


87r2 j 


1/A2 S 






n 


87r2[ 


2 



(2scoths - 2) + 



/ s' 



V sinh^ , 



- 1 ) +£Vs2coths 



cr 



2/72 



(138) 



Here we used relation (H6f) and the relations 



-I3x 



X 



P 



coth x]dx = + - ] - In - - - ] Re/5>0 



P 



(139) 



. .^dr = 2i-^r(/i) 
smh T 



2C(/i-l, f 
/i > 2 



(140) 



All the other functions , F2 can similarly be found 
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APPENDIX B 



In this Appendix we analyze the next-to-leading order in l/Nf. expansion in the 
(3 + 1) -dimensional NJL model (the analysis of the (2 + l)-dimensional NJL 
model is similar). Our main goal is to show that the propagator of the neutral 
NG boson vr has a (3+1) -dimensional form in this approximation and that (unlike 
the (1 + 1) -dimensional Gross-Neveu model [|13|), the l/N^ expansion is reliable 
in this model. 



For an excellent review of the 1/A^c expansion see Ref. For our purposes, it is 
sufficient to know that the perturbative expansion is given by Feynman diagrams 
with the vertices and the propagators of fermions and composite particles a and vr 
calculated in leading order in l/Nc- In the leading order, the fermion propagator 
is given in Eqs. (|13]) and ([T^) (with m replaced by mdyn)- As follows from Eq. 
(P^), the bare Yukawa coupling of fermions with a and vr is Qy"* = 1 in this 



approximation. The inverse propagators of a and vr are [IG]: 



D-p\x) 



A' 



= 1 T 



6\x) + i tT[S{x, 0)TpS{0, x)Tp] 



(141) 



^7^. Here S{x,0) is the fermion propagator 



where p = (a, vr) and To- = 
(|1^) with the mass mdyn = a defined from the gap equation (H). Actually, for 
our purposes, we need to know the form of the propagator of vr at small momenta 
only. We find from Eqs. (|63|) and (|6^): 



Nr 



where 



/l = 


m 


/2 = 


[4- 



+ 1 



1 

&^ 3_ 



(142) 

(143) 
(144) 



The crucial point for us is that, because of the dynamical mass mdyn, the fermion 
propagator (despite its (l + l)-dimensional form) is soft in the infrared region and 
that the propagator of tt has a (3 + l)-dimensional form in the infrared region 
(as follows from Eqs. (^) and (|6^) , the propagator of a has of course also a 
(3 + l)-dimensional form). 

Let us begin the analysis by considering the next-to-leading order corrections in 
the effective potential. The diagram which contributes to the effective potential in 
this order is shown in Fig. la. Because of the structure of the propagators pointed 
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out above, there are no infrared divergences in this contribution to the potential. 
(Note that this is in contrast to the Gross-Neveu modeh the contribution of 
this diagram is logarithmically divergent in the infrared region in that model, i.e. 
the 1/Nc expansion is unreliable in that case). Therefore the diagram in Fig. la 
leads to a finite, 0(1), correction to the potential V (we recall that the leading 
contribution in V is of order Nc). As a result, at sufficiently large values of 
the gap equation in this model admits a non-trivial solution ex 7^ in next-to- 
leading order in l/Nc-, i.e. there is spontaneous chiral symmetry breaking in this 
approximation. 

Let us now consider the next-to-leading order corrections to the propagator of the 
NG mode vr. First of all, note that in a constant magnetic field, the propagator 
of a neutral local field (p{x), D^{x,y), is translation invariant, i.e. it depends 
on [x — y). This immediately follows from the fact that the operators of space 
translations in Eqs. ( p8| ) and (|29| ) take the canonical form for neutral fields 
(the operator of time translations is both for neutral and charged fields in 
a constant magnetic field). The diagram contributing to the propagators of the 
NG mode in this order is shown in Fig. lb. Because of the dynamical mass 
mdyn in the fermion propagator, this contribution is analytic at = 0. Since at 
large the gap equation has a non-trivial solution in this approximation, there 
is no contribution of 0{k^) ~ const, in the inverse propagator of tt. Therefore 
the first term in the momentum expansion in its inverse propagator has the form 
Ci(/i;Q— ^3)— C2k^, where Ci and C2 are functions of a£ and Ai and the propagator 
takes the following form in this approximation: 



(compare with Eq. (|142|) ). 

For the same reasons, there are also no infrared divergences in the fermion propa- 
gator (see Fig.lc) nor in the Yukawa vertices (see Fig. Id) in this order. Therefore, 
at sufficiently large values of Nc, the results retain essentially the same as in lead- 
ing order in l/Nc- 

We believe that there should be no essential obstacles to extend this analysis for 
all orders in l/Nc- 



/i(A£, a£) [(1 - ^Oi(A^, a£)) (fco' - kl) 




(145) 
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APPENDIX C 



In this Appendix we analyze integral equations ( |106| ) and ( |114]) for the function 
A{p). We will show that these equations lead to the expression ( 111| ) for m^y^. 

The integral equations were analyzed by using both analytical and numerical 
methods []. The numerical plot of the dynamical mass as a function of a is shown 
in Fig. 2. Below we consider the analytical analysis of these equations. 

We begin by considering equation (|106|) in the Feynman gauge. After the angular 
integration, this equation becomes 



^ip') = ^J^^K{p^e) (146) 
27r J + miy^ 



with the kernel 



K{p\ e) = I ^^exp(-.fV2) ^^^^^ 



To study Eq.( |146| ) analytically it is convenient to break the momentum inte- 
gration into two regions and expand the kernel appropriately for each region 
(compare with Ref . p] ) : 



a 



271 



deA{e) Jdzexp{-zf/2) 



k"^ + m; 



2 



p^ + z 



+ 



dk^A{k^) Jdzexp{~zf/2) 



k^ + z 



(148) 



Introducing dimensionless variables x = p'^i'^/2, y = k'^i'^/2, a = m'^y^i'^/2, we 
rewrite Eq. ( |148|) in the form 



a 



A(x) = — 
^ ' 2-K 



dyA{y) , r dyA{y)g{y) 



y + a 



y + 



(149) 



where 



dze ^ 
z + X 



-e'Eii-x) 



(150) 



^We thank Anthony Hams and Manuel Reenders for their help in numerical solving these 
equations. 
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and Ei{x) is the integral exponential function. The solutions of integral equation 
( p.49|) satisfy the second order differential equation 



A" - - ^g'^ 
g' 27r X + 



0, 



(151) 



where (/) denotes derivative with respect to x. The boundary conditions to this 
equation follow from the integral equation (|149|) : 



7 

(A 



x=Q 

_g_A^^ 
9' ' 



0. 



For the derivatives of the function g{x) we have relations 



9 = + 9[ 

X 



9 = ^ + 9{x), 



X^ X 



and g{x) has asymptotics 



(152) 
(153) 



(154) 



g(x) ~ In , X — 0, 

X 

112 

g{x) ~ ^ + ^5 x-^oo. 

X x^ x-^ 



(155) 



Using Eqs.( p.54] ),( |155| ) we find that Eq. ( |151| ) has two independent solutions which 
behave as A{x) ~ const, A{x) ~ Inl/x and A{x) ~ const, A{x) ~ ^, near x = 
and X = oo, respectively. The infrared boundary condition (BC) ( |152D leaves 
only the solution with regular behaviour, A{x) ~ const, while the ultraviolet 
BC gives an equation to determine a = a{a). To find analytically a{a) we will 
solve approximate equations in regions x << 1 and x >> 1 and then match the 
solutions at the point x = 1. This provides insight into the critical behaviour of 
the solution at a 0. A numerical study of the full Eq.( p.46| ) reveals the same 
approach to criticality (see the plot in Fig. 2). 

In the region x << 1 Eq. (|151|) is reduced to a hypergeometric type equation: 

A" + -A' + ^ =0. (156) 

X 27rx(x + a^j 

The regular at x = solution has the form 

X / ry 

A,{x) = CiF{tiy, -zu, 1; --), z/ = J—, (157) 

V ZTT 
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where F is a hypergeometric function P|. In the region x >> 1 Eq. (|151|) takes 
the form 



a 



A 



X 2tt x'^{x + a'^) 



0. 



158) 



The solution satisfying ultraviolet BC (|153|) is 



A^ix) = C2-F{^-^, 2; 

X 2 2 X 



2a 



1. 



(159) 



Equating now logarithmic derivatives of Ai and ^2 at x = 1 we arrive at the 
equation determing the quantity a(a): 



d_ 

dx 



In 



xF{iu, -iv; 1; -^) 

2 ' 2 ' .T 



0. 



x=l 



(160) 



Note that up to now we have not made any assumptions on the value of the 
parameter a. Let us seek now for solutions of Eq.( |160| ) with a << 1 (which 
corresponds to the assumption of the LLL dominance). Then the hypergeometric 
function in denominator of Eq.( p.6(J| ) can be replaced by 1 and we are left with 
equation 



1 



1 



1 



— z/'F(l + iu,l- iv- 2; -) + F{iv, -iu; 1; 



0, 



(161) 



where we used the formula for differentiating the hypergeometric function P| 

(162) 



d „, , ^ ab 

-rF(a, b, c; z) = —Fia + 1, 6 + 1, c + 1; z). 

dz c 



Now, because of a << 1, we can use the formula of asymptotic behavior of 
hypergeometric function at large values of its argument z p[: 



^, , , r(c)r(6-a) 

F{a,b,c; z) ~ t[-z] 



V{c)T{a-b)^ 



r(6)r(c-a)' ' r(a)r(c-6) 

Then Eq.( |161| ) is reduced to the following one: 



(163) 



cos 



t/ln — + argS(z/) 



l + ivT{l + 2iu) 
2 T^il + iv) 



and we get 



dyn = 2|e5| exp 



TT 



(2r2 + l)/2 -argS(z/) 



(164) 



(165) 
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where n is zero or positive integer. The argE(z/) can be rewritten as 

argS(z/) = arctani/ + argr(l + 2iv) — 2argr(l + iv) (166) 
and in the hmit z/ Eq. (|160|) takes the form 



m 



dyn 



2|e-B|eexp 



TT 



-(2n+l) 



2|ei?|eexp 



TT 



(167) 



(the second factor e here is e ~ 2.718 and not the electric charge!). The stable 
vacuum corresponds to the largest value of m\-^^ with n = 0. 



Let us now turn to equation (|114|) in an arbitrary covariant gauge. The function 
g{x) is now replaced by 



^- ^ = g{x) + -\xg\x). 

z + X 2 



168) 



As it is easy to verify, this does not change equation ( |156| ) in the region a; << 1. 
At X » 1 we have 



g'{x) ~ 



2-X _^1-A_ 



2x2 



22- A-6 



l-A 



g'[x) 



x2- A-4 



l-A • 



(169) 



Therefore in any gauge, except A = 2, the different ional equation at x >> 1 takes 
the form 



A 



X Arc x^lx + 







;i70) 



with asymptotic solution A(x) ~ ^. In the gauge A = 2, instead of Eq. ( |170D , we 
have 



A" + -A' + 3 

X 7rx'^(x + a^j 



0, 



(171) 



which gives more rapidly decreasing behavior A{x) ~ ^ when x —>■ oo. Repeating 
the previous analysis, we are led to expression ( |111| ) for mdyn- 



49 



References 

[1] J. Schwinger, Phys. Rev. 82 (1951) 664. 

[2] V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Phys. Rev. Lett. 73 
(1994) 3499; Phys. Rev. D 52 (1995) 4718. 

[3] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345. 

[4] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 108 (1957) 1175. 

[5] I. S. Gradshtein and I. M. Ryzhik, Table of Integrals, Series and Products 
(Academic Press, Orlando, 1980). 

[6] V. A. Miransky, Dynamical Symmetry Breaking in Quantum Field Theories 
(World Scientific Co., Singapore, 1993). 

[7] B. Simon, Ann. Phys. 97 (1976) 279. 

[8] A. I. Akhiezer and V. B. Berestetsky, Quantum Electrodynamics (Inter- 
science, NY, 1965). 

[9] A. CappeUi, C. A. Trugenberger, and G. R. Zemba, Nucl. Phys. B396 (1993) 
465. 

[10] A. Chodos, I. Everding, and D. A. Owen, Phys. Rev. D42 (1990) 2881. 

[11] T. Appelquist, M. Bowick, D. Karabali, and L. C. R. Wijewardhana, Phys. 
Rev. D33 (1986) 3704. 

[12] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17 (1966) 1133; S. Coleman, 
Commun. Math. Phys. 31 (1973) 259. 

[13] D. Gross and A. Neveu, Phys. Rev. DIO (1974) 3235; E. Witten, Nucl. Phys. 
B145 (1978) 110. 

[14] J. Zak, Phys. Rev. 134 (1964) A1602; J. E. Avron, I. W. Herbst, and B. 
Simon, Ann. Phys. 114 (1978) 431. 

[15] S. P. Klevansky, Rev. Mod. Phys. 64 (1992) 649. 

[16] V. P. Gusynin and V. A. Miransky, Mod. Phys. Lett. A6 (1991) 2443; Sov. 
Phys. JETP 74 (1992) 216; V. A. Miransky, Int. J. Mod. Phys. A8 (1993) 
135. 

[17] D. Caldi, A. Chodos, K. Everding, D. Owen, and S. Vafaeisefat, Phys. Rev. 
D39 (1989) 1432. 



50 



[18] P. I. Fomin, V. P. Gusynin, V. A. Miransky, and Yu. A. Sitenko, Riv. del 
Nuovo Cimento 6, N5 (1983) 1. 

[19] A. M. Perelomov and V. S. Popov. Theor. Math. Phys. 4 (1970) 664. 

[20] W. Dittrich and M. Renter, Effective Lagrangians in Quantum Electrody- 
namics ( Springer- Verlag, Berlin, 1985). 

[21] G. Calucci and R. Ragazzon, J. Phys. A27 (1994) 2161. 

[22] J. Schwinger, Phys. Rev. 125 (1962) 397; J.H. Lowenstein and J.A. Swieca, 
Ann. Phys. (N.Y.) 68 (1971) 172. 

[23] R. Parwani, Phys. Lett. B358 (1995) 101. 

[24] D. Cangemi, E. D'Hoker, and G. Dunne, Phys. Rev. D51 (1995) R2513; 
Phys. Rev. D52 (1995) R3163. 

[25] R. Jackiw, Phys. Rev. D29 (1984) 2375; 1. V. Krivc and A. S. Rozhavsky, 
Sov. Phys. Usp. 30 (1987) 370; A. Kovner and B. Rosenstein, Phys. Rev. 
B42 (1990) 4748; G. W. Semenoff and L. C. R. Wijewardhana, Phys. Rev. 
D45 (1992) 1342; N. Dorey and N. E. Mavromatos, Nucl. Phys. B368 (1992) 
614. 

[26] Y. Hosotani, Phys. Lett. B319 (1993) 332. 

[27] T. Vachaspati, Phys. Lett. B265 (1991) 258; K. Enqvist and P. Olcscn, Phys. 
Lett. B319 (1993) 178; J. M. Cornwall, in Unified Symmetry: In the Small 
and Large, eds. B. N. Kursunoglu et al. (Plenum Press, New York, 1995). 

[28] N. K. Nielsen and P. Olesen, Nucl. Phys. B144 (1978) 376. 

[29] P. Salabura et al, Phys. Lett. B245 (1990) 153; I. Koeing et al, Z. Phys. 
A346 (1993) 153. 

[30] L. S. Celenza, V. K. Mishra, C. M. Shakin, and K. F. Lin, Phys. Rev. Lett. 
57 (1986) 55; D. G. Caldi and A. Chodos, Phys. Rev. D36 (1987) 2876; Y. 
J. Ng and Y. Kikuchi, Phys. Rev. D36 (1987) 2880. 

[31] D. G. Caldi and S. Vafaeisefat, Phys. Lett. B287 (1992) 185. 

[32] S. Coleman, Aspects of Symmetry (Cambridge University Press, Cambridge, 
1985). 



51 



FIGURE CAPTIONS 



1. Fig. 1. Diagrams in next-to-leading order in l/N^- A solid line denotes the 
fermion propagator and a dashed line denotes the propagators of a and tt 
in leading order in 

2. Fig. 2. A numerical fit of the dynamical mass as a function of the coupling 
constant a: m^y^j \/1eB — exp[— |y'^a -|- 6]. The fitting parameters are: 
a = 1.000059, h = -0.283847. 
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a) 



b) 





c) 



Figure 1: Diagrams in next-to- leading order in l/Nc- A solid line denotes the 
fermion propagator and a dashed line denotes the propagators of cr and tt in 

leading order in l/N^. 
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Figure 2: A numerical fit of the dynamical mass as a function of the coupling 
constant a: m^y^j \/1eB — exp[— |y'^a + b]. The fitting parameters are: 

a = 1.000059, -0.283847. 
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